
clc;
clear;
%% 条件准备
% 首先设置涡旋nc文件路径
AEfile = 'H:\ftp\eddy\anticyc1993_2021_2cm.nc';
% 设置范围
extent = [0 25; 108 122];
% 设置时间
year = [1993 2021];



% 可以采取两种方法初始化类
% 1.通过构造函数,然后调用set方法挨个添加属性,set方法可以自行添加
EDC = EddyDatasetClass(AEfile);
EDC.extentSet(extent);
EDC.yearSet(year);

% 2.通过构造方法的varargin不定输入进行初始化，可以自行设定
% EDC = EddyDatasetClass(file, 'extent',LNW_extent,'year',[1993 2021]);


%% 涡旋文件读取与处理

% 从nc文件中读取到Dict中
EDC.eddyFile2Dict();

% 将Dict数据转化成轨迹数据
EDC.Dict2traj();

% % 将Dict数据转化为eddyArray涡旋数组
% eddyArray1 = EddyDatasetClass.eddyArrayAddByDict([], EDC.Dict, -1);
% 
% % 将Traj数据转化为eddyArray涡旋数组
% eddyArray2 = EddyDatasetClass.eddyArrayAddByTraj([], EDC.traj, -1);


%% 季节分割
% spilt_traj = EddyDatasetClass.trajSplitBySeason(EDC.traj,[3 4 5;6 7 8;9 10 11;12 1 2]);

save('AE_EDC_2.mat','EDC');
